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Routers, switches, and repeaters are essential components of modern information-processing systems. 
Similar devices will be needed in future superconducting quantum computers. In this work we investigate 
experimentally the time evolution of Autler-Townes splitting in a superconducting phase qubit under the 
application of a control tone resonantly coupled to the second transition. A three-level model that includes 
independently determined parameters for relaxation and dephasing gives excellent agreement with the 
experiment. The results demonstrate that the qubit can be used as a ON/OFF switch with 100 ns operating 
time-scale for the reflection/transmission of photons coming from an applied probe microwave tone. The 
ON state is realized when the control tone is sufficiently strong to generate an Autler-Townes doublet, 
suppressing the absorption of the probe tone photons and resulting in a maximum of transmission. 
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In the optical regime, a variety of effects that can be employed for switching have been demonstrated, for 
example photon blockade by single atoms \ parametric instabilities in nonlinear optical wave-mixing^, her- 
alded single-photon absorption by one trapped ion^, single-atom dynamic control of light in microresonator 
photon turnstiles^, and population inversion in single dye molecules^. Optical switching using the phenomenon of 
electromagnetically induced transparency^ has been realized in gases of magneto -optically trapped ^^Rb^; recently, 
it has become possible to use for this purpose only a single ^^Rb atom in a high-fmesse cavity^. With improved 
dissipation and decoherence, it could be possible to use electromagnetic-induced transparency effects to build 
similar devices with superconducting quantum circuits^. Other applications employing multilevel systems - for 
example state preparation, emulation of large quantum-number spins^°, single-shot fast quantum nondemolition 
readout techniques^\ implementation of controlled Z-n gates in certain quantum computing architectures^^, 
single-qubit microwave ampliflers^^ etc. - can be envisioned. 

In this paper we demonstrate the dynamical operation of a phase qubit as a microwave photon-absorbing 
device that can be switched ON or OFF in approximately 100 ns by using an external control field. The function- 
ing of our device is based on an effect closely related to electromagnetically induced transparency: under the 
application of a relatively intense continuous field (the control field) to the second transition of the qubit, the 
spectral line of the first transition splits into a doublet. This is called Autler-Townes effect^^, and it has been 
recently observed in superconducting quantum systems as welP^"^^. Here we operate our device in time-domain, 
which allows us to observe transient effects that occur until a steady state is reached. We observe how the doublet 
forms when the control field is suddenly applied, and how the spectral line of the first transition appears back 
from the merging of the Autler-Townes peaks when the coupling field is switched off. 

Besides the fundamental importance for understanding the dynamics of the Autler-Townes effect, our experi- 
ment has potential applications in the field of microwave photonics, for example for fast quantum switches that 
can be integrated in superconducting quantum-processor architectures. This can be achieved by realizing that the 
"control" microwave tone applied to the qubit determines the absorption rate for a "probe" microwave tone. 
Indeed, when the control field is OFF, the third level does not play any role and the device absorbs the incident 
probe field radiation resonant to the lowest transition, and its first excited level is populated. If such a two-level 
system is embedded into a transmission line, it has been predicted theoretically in quantum optics^^ and demon- 
strated experimentally^^ that the qubit will reflect back the incoming photons by resonant fluorescence, provided 
that the intensity of the probe field is not too high (this saturates the qubit). When the control field is ON, the 
absorption from the qubit is suppressed almost completely - the system becomes "transparent" to the incident 
probe radiation and the transmission is maximal. We calculate theoretically the expected reflection coefficient 
both in the transient and in the steady-state regimes, and we also analyse the performance of the switch at different 
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strengths of the control and probe tones. Using our theoretical 
model, we are able to present a set of general guidelines for the design 
of superconducting switches based on the Autler-Townes effect - we 
show that high anharmonicity, low probe powers, and low dephasing 
rates lead to an effective suppression of the reflection. 

Results 

Our phase qubit^°'^^ consists of an r.f. SQUID with loop inductance L, 
junction capacitance C, and Josephson energy £j, which is read out by 
another d.c. SQUID fabricated on-chip in the proximity of the qubit 
and coupled inductively to it (see Fig. 1). The qubit can be biased by 
an externally applied magnetic flux and can be coupled to other 
circuit elements such as transmission lines and resonators. The 
device forms a multilevel quantum system with the first three levels 
denoted by |0), and |2). These energy levels can be addressed 
using microwave fields closely resonant to the first |0) ^ 1 1) and the 
second |1) ^ |2) transitions (the |0) ^ |2) transition is almost for- 
bidden). 

A detailed description of the functioning of the device has been 
given elsewhere^^. The Hamiltonian of the system is 

H= \- EjC0S\2n--], (1) 

2C 2L ^ \ ^oj 

where Q is the charge on the capacitor formed by the junction, ^ is 
the flux variable, % = h/2e = 2.067X 10"'^ Wb is the magnetic flux 
quantum, and Ogxt = ^dc+^rf (0 is the total magnetic-flux com- 
ponent of the externally applied d.c. and r.f fields. In this system, the 
r.f tones applied have angular frequencies Wp (probe tone) and Wc 
(control tone) close to the two resonant frequencies coio of |0) ^ 1 1) 
transition and (021 of |1) ^ |2) transition, respectively (see Fig. 2). 
The first two transitions were determined in independent single-tone 
and two -tone spectroscopy measurements^^, and they had the values 
Wio = 2n X 8.135 GHz and CO21 = 2n X 7.975 GHz. The probe and 
control fields are supplied to the device through an on-chip coplanar 
waveguide transmission line. The Rabi frequencies of these fields are 
denoted by Qp and respectively, and they are directly proportional 
to the corresponding probe and control field amplitudes. 

Experimental data. Here we demonstrate that this three-level 
superconducting device can be operated dynamically, in time 
domain. In Fig. 2 we present the pulse sequence applied in order 
to use the system as a switch based on the Autler-Townes effect. 
In Fig. 3 we show the results of the experiment. The ON state is 
characterized by the suppression of excitations to the state |1) 
normally caused by absorption of photons from the probe field. It 
is produced by a large value of the control field amplitude (effective 
Rabi frequency Q^)- The state OFF corresponds to a low enough 
value of the control field so that excitations to the state |1) are 
allowed. Although we do not measure the transmission (reflection) 
of photons directly, the data demonstrates unambiguously that the 




Figure 1 | Image of the Josephson phase qubit sample used in the 
experiment. An Ml MO J M junction with an area of ~ 14 fim^ was 
patterned on a sapphire substrate by standard optical lithography. SiN;^ 
was used as an insulator between different metallic layers. 




delay 

Figure 2 | Schematic of the pulse sequence for the dynamic Autler- 
Townes effect. The qubit is probed by a continuous microwave field 
{Wpy Qp), with the frequency Wp being swept around the first transition 
frequency Wiq. Each measurement sequence starts at t = 0 ns with the 
system in the stationary state as given by the application of the probe field. 
At t = 150 ns we apply the coupling pulse (co^, Qc)> exactly resonant with 
the second transition, €0^= CO21. 

device can be used to modulate the transfer of photons by controlling 
the absorption rate of the qubit. The measured state occupation 
numbers show clearly the dynamical process of formation of the 
Autler-Townes doublet, how it reaches the stationary state, and 
how the two peaks finally collapse into a single spectroscopic 
signal at the first transition frequency when the coupling field is 
switched off. In both the on-set and the switch-off stages of the 
effect, the time- scale for the system to reach the steady state is of 
the order of 100 ns, in agreement with the decoherence times in our 
system (see the theoretical model below). 

Numerical simulation. In order to simulate the experimental results 
including the effects of decoherence, we use the standard Liouville 
superoperators in the Markov approximation. For our system, 
described by the Hamiltonian equation (1) and placed in a 
reservoir of electromagnetic modes at zero temperature (see^^ for a 
detailed derivation), the Markovian master equation for the density 
matrix p in the rotating frame reads 

p=-^[H,p]+Ael[p]+^depW. (2) 

n 

The effective Hamiltonian H is obtained by truncating equation 
(1) to the lowest three levels and performing a rotating wave 
approximation. It can be put in the form 

H=^ [ApGu + (Ap + A,)(722] + 

+ ^[(Qpe^'^ + 0.69Q,e-^'^>oi 

+ (Q, + 1.45Qpe^'^^+^'^)(7i2]+h.c.. 

Here we have taken an external r.f magnetic field of the form Oj-KO = 
Op cos{cOpt + (/>) + cos{cDct). This produces an effective coupling 
to the qubit g{t) = gp cos{cOpt + (j)) + gc cos{cOct), where 
gc=- i^c/L) y/l/Chcoo and gp=- {%/ L) ^I/C/zcoq (where Wq 
is the Josephson plasma frequency) resulting in Rabi frequencies 
Qp = 0.69gp and — gc (see Ref 22 for a detailed derivation). 
Here we have used the operators aij = \i){j\ (which generalize the 
Pauli matrices to a many-level system), and we have defined 
the detunings Ap = Wiq — cOp, = CO21 ~ ^rid 3 = Wp — (0^. 
The origin of the qubit-field coupling term [the second line in the 
Hamiltonian equation (3)] can be understood easier if one treats the 
device as a harmonic oscillator interacting with fields; in this 
case, the corresponding terms in the Hamiltonian read 
{h/2) [{I/V2) {gpe^^ +gce-^'')cJoi + {gc+gpe'''^'^)c7u] +h.c.. The 
result of including anharmonic corrections is that from the 
previous expression is replaced by 1.45 and l/v5 by 0.69. 
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Figure 3 | On-set and extinction of the Autler-Townes effect. The figure shows the experimentally determined sum of excited states level-occupancy 
probabilities, P\ i) + P|2>, obtained by applying the pulse sequence from Fig. 2. The arrow indicates the frequency of the first transition, 
coio/271 = 8.135 GHz. 

Fio and r2i, respectively, y^Q, are intrinsic (pure) dephasing rates 
of the states |1) and |2), and k = \JY\^Y2\. The dissipation para- 
meters were either measured directly in independent experiments 
{e.g. by exciting the qubit and measuring the decay time) or extracted 
from spectroscopy data^^'^^, and for this sample they were = 271 X 
1 1 MHz, Tio = 271 X 7 MHz, = 2nx7 MHz, y%=2nxl6 MHz. 

In Fig. 4 we present the result of the numerical simulation using 
equation (2). For this simulation, the amplitude of the probe field 
corresponds to Rabi frequency Qp = 2nX 3.45 MHz; for the control 
field, the corresponding values were — 271 X 60 MHz in the ON 
state, and a much smaller value = 27i X 5 MHz in the OFF state. 
The calibration of the probe and control field amplitudes was done in 
independent measurements of the Rabi frequency between the levels 
|0) and and between the levels |1) and |2) respectively. In the 
absence of the control field, the level |2) is empty while the level |1) 
is populated when the probe field is resonant with the |0) ^ |1) 




8.05 '1 00 

Figure 4 | Simulated occupation probabilities for the sum of the occupation probabilities for the states 1 1) and |2) with averaging over the phases ^. 

This quantity is measured directly in the experiment (see Fig. 3). 



The relaxation part is given by^^ 

>^rel[p] = ^(2o-oipO-io-0-iip-p(7ii) 

, ^21 X (4) 

+ ^(2^12^(721 -a22p-p(722) ^ ^ 

+ K:(e"^^^(7oiP(72l+e'^Vi2p(7io), 

and the dephasing is 

+ ^(2(722PO-22-0-22p-pO-22). 

Here p is the density matrix in a doubly- rotating frame^^, the inter- 
level relaxation rates between 1 1) ^ |0) and |2) ^ 1 1) are denoted as 
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transition. After the control field has been switched on, the original 
high occupation peak centred at coio splits into the Autler-Townes 
doublet. These observations are in good agreement with the experi- 
ment (Fig. 3). The simulation also captures two additional features of 
the experimental data: the asymmetry in the Autler-Townes peaks 
and the appearance of a plateau with non-zero occupation probabil- 
ities even when cjp is far off-resonant from coio in the ON state. Both 
of these effects result from the existence of cross-coupling terms in 
the Hamiltonian equation (3), see the Discussion below. 

In the experiment, the phase difference between the probe and 
the control field cannot be maintained at the same value for each 
realization of the sequence presented in Fig. 2. To model this, we 
calculate the probabilities P|i> + P\2) for several values of (j) randomly 
distributed in the interval [0, 2n] and we average over the results. 
More considerations about the effect of the cross-coupling terms are 
delegated to the Discussion. 

We note that the power levels used in this experiment correspond 
to operation close to the level of single microwave photons. To 
estimate the average number of photons involved, consider a generic 
process (without dissipation) of resonant Rabi oscillations with Rabi 
frequency Qr. The probability for a two-level system to be excited in a 
time T is sin^(QRT/2). Therefore, after a time Tr = ti/Qr the system is 
with certainty in the first excited state, meaning that it has absorbed 
exactly one photon. The corresponding power is ^co/tr. In our case, 
the probe Rabi frequency Qp= 27i X 3.5 MHz results in Xp = 143 ns 
(and a power of fiOJiolip = 38 aW = — 134 dBm), and for the control 
field in the ON state, 0.^= 27i X 60 MHz yields = 8.3 ns (and a 
power of hoJiJ^c — 0.6 fW = —122 dBm). Thus, during the time 
between t = 150 ns and t = 300 ns there were on average about 1 
photon present in the probe beam and about 18 photons in the 
control beam. Moreover, as shown in the next section, increasing 
the input power level saturates the qubit and switching does not 
occur anymore. These low power levels are typical for measurements 
of superconducting qubits, therefore the devices based on the mech- 
anism described in this paper can be seamlessly integrated in any 
circuit QED architecture. 

Discussion 

We start by discussing the effect of the cross -coupling terms and that 
of the relative phases between the coupling and the probe tones. In 
the effective Hamiltonian equation (3), there are oscillating terms 
with factors exp(±i3t), which result from cross couplings of the 
coupling field into the |0) ^ |1) transition and of the probe field 
into the 1 1) ^ |2) transition (see also the detailed derivation in Sec. II 
and Sec. Ill of Ref 22). In the real experiment, it was not possible to 



ensure that the phase between the coupling and the probe fields 
remained the same for every run of the measurement sequence. To 
obtain Fig. 4, we have averaged over which is what occurred 
effectively in the experiment as well. It is however instructive to see 
what would happen for a constant phase, and for the case in which 
there is no cross-coupling at all. 

In Fig. 5(a) we show the result of simulations with the same master 
equation as equation (2) but with the relative phase between the 
probe and the coupling fields fixed at (/> = 0. One can see the appear- 
ance of fringe-like structures in the spectrum. This figure can be 
regarded as a snapshot corresponding to a single run of the pulse 
sequence. The plateau seen in the previous spectra in the region 
where the coupling field is on and the probe is far off-resonant still 
exists; also the asymmetry between the Autler-Townes peaks is pre- 
served. At different (/>'s, these fringes appear in different places; thus 
when averaging over (j) they will be washed away. 

To get even more understanding of the effects of the coupling 
terms, let us analyse what happens if there is no cross -coupling at 
all. This would be the case for a strongly anharmonic system. In this 
case we can approximate the effective Hamiltonian as 

H'=^ [Ap(7ii + (Ap + A,)(722] + ^ [QpGoi+^cC^u] +h.c. (6) 

This Hamiltonian is widely used in quantum optics for studying 
electromagnetically induced transparency and coherent population 
trapping^. 

Fig. 5(b) shows a time-domain simulation similar to that in 
Fig. 5(a) but with the Hamiltonian H'. The timeoscillating terms in 
£i-ei [p] are similarly neglected. Several features can be noticed imme- 
diately. The Autler-Townes peaks are this time symmetric, the occu- 
pation probabilities at cOp far away from Autler-Townes splittings are 
smaller, and we do not have a plateau anymore in the ON state (from 
150 ns to 300 ns). The latter two features can be explained intuitively 
in a simple manner: the coupling field tends to introduce an addi- 
tional population of level |1) when it couples into the |0) ^ |1) 
transition. This effect does not depend on the frequency of the probe 
field, which is why it is appears as a plateau. 

In the remaining part of this section, we will discuss the operation 
of the phase qubit as a microwave quantum switch. The capability of 
distributing entangled microwave photons in superconducting 
quantum networks is an essential ingredient for building future 
quantum processors'^ An immediate application of the effect pre- 
sented above is to realize a quantum switch for controlled routing of 
microwave photons. These devices can be realized by embedding the 
qubit into an open transmission line and using a vector network 
analyser to measure the scattering coefficients of the probe field. 



(a) 



Occupation probabilities P|^^+ P|2>(simulation, (\>=0) 



(b) 



Occupation probabilities P|i>+ P|2>(without cross-coupling) 




400 



8.15 

probe frequency 
(0 /27r(GHz) 8.1 



200 

measurement 
delay (ns) 



8.05 




8.15 

probe frequency 
CO / 271 (GHz) 8.1 



200 

measurement 
delay (ns) 



8.05 



Figure 5 | Simulated + P|2> with ^ = 0, and (a) with cross-coupling, (b) without cross-coupling. 
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yielding the corresponding reflection and transmission coefficients. 
In the stationary regime, these measurements have been done for flux 
qubits^ and transmons^^. Here we discuss the general dynamic char- 
acteristics of these devices, with exemplification by using the para- 
meters corresponding to our qubit. 

For qubits embedded into one-dimensional transmission lines, the 
reflection coefficient r of the probe tone is proportional to the off- 
diagonal element pio = Tr[p(7io]^: 



(7) 



where rj is the dimensionless coupling efficiency to the transmission 
line^^. In the case of our experiment, the switching occurs between 

(st) (si) 

two Steady states, and Pqff' which are reached after —100 ns of 
transients. As a figure of merit for switches, the OFF/ON power ratio 
between the reflected powers {R = \ r/rj\^) in the OFF state and that in 
the ON state i^oFF/^oN must be as large as possible. In the simula- 
tions below we will use parameters corresponding to our experiment; 
however, we stress that the theoretical results are generic and can be 
applied to any type of superconducting quantum multi-level system. 
Below we will analyse systematically the influence on the quantum 
switch performance of the following effects: cross-coupling (anhar- 
monicity), saturation, and dephasing. 

We first calculate numerically the off- diagonal density matrix 
element p^Q and the reflection coefficient obtained from equation 
(7) in the steady state. For these calculations we take both the probe 
and control fields resonant with the respective transitions (Ap = 0 
and Ac = 0). This corresponds to the experimental data along the 
arrow from Fig. 3. In Fig. 6(a) and (b) we present the normalized 
reflected power R = \r/rj\^ simulated with cross-coupling (using the 
Hamiltonian H), and respectively without cross-coupling (using the 
Hamiltonian H'). The effect of the cross-coupling terms become 
apparent especially at high values of Q^; this field couples into 
the |0) — |1) transition, the reflected power is not turned down 
completely. We have checked numerically that this is due to the real 
part of PiQ acquiring a non-zero value, while the imaginary parts of 
p^^Q with and without cross -coupling do not differ much. For 

example, without cross-coupling, Rej^pJo^'j =0, Im j^pJo^'j ~ — 0.19 
in the OFF state, and Re[p[o^] =0, Im[p[o^] ^ -0.02 in the ON 
state, which give a power ratio Rqff/Ron of 90; whereas with 
cross-coupling. Re [pJo^] ^ -0.01, Im^pJ?] ^ -0.19 in the OFF 



state, and Re[pJo^] ^ -0.1, Im^pJo^j ^ -0.04 in the ON state, 

which give a ratio Rqw^Ron of only 3. 

This means that, in the ON state, the normalized power reflection 
R with cross-coupling is always larger than that without cross-coup- 
ling, and the additional power reflection due to the cross -coupling 
reduces the ON/OFF ratio, Rovv^RoNy of the switch. We conclude 
that, in order to have an efficient switch, it is advantageous to 
increase the anharmonicity of the qubit, which would reduce the 
cross -couplings [for example in equation (3), if 3 is large, the 
cross-coupling terms will oscillate fast and average to zero]. Some 
types of qubits, for example the flux qubit, do have a high degree of 
anharmonicity by design; for phase qubits there exist as well propo- 
sals for increasing the anharmonicity^^. 

If the probe power is too large, from Fig. 6 we can see that, in the 
OFF state, the power reflection R is significantly reduced. This is 
because in the OFF state the system is essentially a two -level system 
exposed only to the probe field radiation. In this case, if Qp is much 
larger than yio, the stationary state of the system in the {|0), |1)} 
subspace is maximally mixed, and the reflection coefficient drops 
to zero. This shows that increasing the power of the probe field would 
not improve the characteristics of a switch. 

In the absence of cross-couplings and for small enough values of 

the probe field, (Qp/y^QY^l and (Qp/Qc)^^l> analytical results 

(st) 

for for the off-diagonal element p^g the steady-state are avail- 
able^^ 



2A, 



- 4 a; + + 710720 + {y,^ + y^o) 



(8) 



where y^g ^ ^^lo + Yio = 2nxl4 MHz and ^ ^21 + 720 = 2nx27 
MHz. In Fig. 7, the real and imaginary parts of the reflection coef- 

(st) 

ficient obtained with p\q of equation (8) are plotted for different 
values of Q^- 

From equation (8) we can also obtain an explicit expression for the 
reflection coefficient r in the steady state at probe-field resonance, 

^10720 



-ri 



^c+y 10720 



(9) 



Using equation (9) we find that in the OFF steady state the reflection 
coefficient is Vq^^^ rj ~ 0.47. Ideally, the refiection coefficient could 
approach the unity, but for the parameters corresponding to our 
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qubit this does not happen due to the relatively large dephasing rate, 
YiQ ~ 2rio- Now, when the coupling field is ON and the system has 

reached the stationary state Pq^, equation (9) predicts that there is 

almost no reflection (^Vq^ < 0.05^ of the probe field (green circle 

in Fig. 6). These results are in reasonably good agreement with the 
numerical simulations in time domain (see also Fig. 8), even though 
the approximation (Qp/Qc) ^1 is not valid stricto sensu. For 
example, equation (9) predicts that the ON/OFF ratio i^oFp/^ON is 
of the order of 100, in agreement with the numerical result obtained 
before. To obtain a higher reflection in the OFF state, equation (9) 
shows that one needs to have a lower pure dephasing rate and to 
reduce the value of the coupling fleld in the OFF state. If the coupling 
fleld in the OFF state is well below the value of the decoherence rates, 

we have from equation (9) that i^oFF / ^on^ [l + ^^c(on) / (^10720)] • 
Experimentally, lowering the OFF field to zero is limited by the 
leakage of the mixer used to pulse the coupling field, while the reduc- 
tion of pure dephasing is a general problem in the fleld of super- 
conducting qubits. 




measurement delay (ns) measurement delay (ns) 



Figure 8 | (a) Real part, (b) imaginary part, and (c) Smith chart of the 
normalized reflection coefficient in time domain, (d) The corresponding 
(normalized) power reflection. The control field is OFF before the 150 ns 
delay time (red), then switched ON between 150 and 300 ns (green), and 
finally switched OFF again (blue). The dashed arrows in (c) indicate the 
direction of time- evolution. The simulations correspond to the 
experimental data along the arrow from Fig. 3. 



In Fig. 8 we present the results of time-domain numerical simula- 
tions of the reflection coefficients for a switching event correspond- 
ing to the one occurring in the experiment presented in Fig. 3. We see 
that, as in the case of the population number, the transients tend to 
stabilize after times of the order of 100 ns. 

To conclude, our results demonstrate that a superconducting 
three-level system with large enough anharmonicity and low dephas- 
ing can be used for example as the node of a quantum network, 
redirecting single photons along chosen paths (reflecting them in 
the OFF state and allowing them to pass in the ON state). 

Methods 

To probe the qubit we use a weak continuous probe tone. The Hamiltonian equation 
(1) is dynamically manipulated by changing the amplitude of the control field, 
resonantly coupling the two higher levels 1 1) and |2). This field is turned on at 150 ns 
and turned off at 300 ns. For measuring the occupation probabilities, we use the fact 
that the three levels reside in a metastable welP°'^\ The potential barrier that creates 
this well is then lowered by a measurement pulse applied after a measurement delay 
time (see Fig. 2), such that the states above the ground state are allowed to tunnel out 
of the well with a certain tunnelling probability. The resulting variation in the mag- 
netic flux piercing the qubit loop is detected by a nearby d.c. SQUID, which is 
switched in the running state^*^. Therefore from the switching statistics we can cal- 
culate the total state occupation probability of the two excited levels P| 1) + P\2f^'^^ that 
have tunnelled out. The numerical simulations are done by solving the master 
equation with the fourth-order Runge-Kutta method^^. 
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